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ABSTRACT 

We describe a uniform all-sky survey of bright blazars, selected primarily by 
their flat radio spectra, that is designed to provide a large catalog of likely 7- 
ray AGN. The defined sample has 1625 targets with radio and X-ray properties 
similar to those of the EGRET blazars, spread uniformly across the |6| > 10° sky. 
We also report progress toward optical characterization of the sample; of objects 
with known R < 23, 85% have been classified and 81% have measured redshifts. 
One goal of this program is to focus attention on the most interesting (e.g., high 
redshift, high luminosity, ...) sources for intensive multiwavelength study during 
the observations by the Large Area Telescope (LAT) on GLAST. 

Subject headings: BL Lacertae objects: general — galaxies: active — quasars: 
general — surveys 



1. Introduction 



It is well-known ( jHartman et al.l Il999l ; iMattox et al.l l200ll ) that many of the high 



latitude EGRET sources are associated with the bright, flat radio spectrum AGN known 
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as blazars. Sowards-Emmerd et al. (2003; SRM03) quantified such associations, develop- 
ing a combined figure of merit (FoM), which measured the likelihood tha t an individual 
radio /X-ray source near the large (~0.7°) Third EGRET Catalog (3EG, lHartman et al. 



19991 ) uncertainty regions is the 7-ray counterpart. They also noted that there are many 
radio-loud blazars with very similar properties not obviously associated with a 3EG source. 
A likely explanation is that blazars are very vari able at high energy, with duty cycles fo r 
the bright, flaring state as small as a few percent (lHartman et al.lll993l ; iKniffen et al.lll993l ). 
During the limited (typically two weeks per pointin g direction) 3EG exposu r e man y of these 
sources may have been in quiescence. Accordingly, ISowards-Emmerd et al.l (120051 ) extended 
the SRM03 analysis by selecting "3EG-like" blazars, i.e., sources whose radio flux density 
and spectrum (and X-ray flux) were very like those of the 3EG blazars but which happened 
not to lie within a 3EG test statistic (TS) uncertainty region. The positions of these sources 
showed a clear excess of 7-ray photons over background and these sources are likely to show 
7-ray high states during future missions. 

The Large Area Telescope (LAT) on GLAST will provide an improvement of several 
orders of magnitude over EGRET/ CGRO with an increased sensitivity in the 50 MeV — 300 
GeV range and a wide (>2.5 sr) field of view. The LAT should detect many thousands of 
sources during the 5- to 10-year mission, with a large fraction of the high-latitude sources 
being blazars. The early mission will be devoted to a sky survey, covering the entire sky 
at good sensitivity every three hours. This will greatly enhance the likelihood of detecting 
transient and variable sources, such as blazars. While several larg e samples of blazars have 
been compiled rece ntly (see especially th e ASDC blazar catalog, iMassaro et al.l 120071 . and 
the ROXA catalog, iTurriziani et al.l 120071 ). there is a surprisingly incomplete knowledge of 
the radio-bright, flat-spectrum population, which is most clearly associated with the GeV 
7-ray sources. We seek to rectify this by defining CGRaBS, the Candidate Gamma- Ray 
Blazar Survey, a large sample of EGRET-like blazars selected across the extragalactic sky. 
By obtaining optical classifications and redshifts for a large fraction of these sources, we 
plan to enable prompt, intensive follow-up of the most interesting (e.g., high redshift, high 
luminosity, peculiar spectrum) sources that are detected in the LAT sky survey data. Fur- 
thermore, identification of a substantial fraction of the LAT sources with blazars will allow 
us to focus on the non-blazar remainder, potentially isolating new classes of cosmic 7-ray 
emitters. 
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2. Sample Selection 

For any FoM-type count erpart selection, it is important to have uniform parent pop- 



ulations. iHealey et al.l (120071 ) have re cently developed s uch a catalog, CRATES, which ex- 



tended results of the CLASS survey jMvers et al.l l2003h to obtain 8.4 GHz observations of 



all |6| > 10° objects brighter than 65 mJy at 4.8 GHz with spectral indices a > —0.5 (where 



S oc y a ). To estimat e the radio spectr al index of the core, we use the NVSS (ICondon et al. 



19981 ) and SUMSS (IBock et al.l 119991 ) lower-frequency surveys. The result is a sample of 
over 11,000 flat-spectrum radio sources with interferometric measurements at ~1 GHz and 
8.4 GHz (with FWHM beam sizes ~40" and ~0.25" respectively), giving precise positions, 
spectral indices, and morphologies for the compact components. The CRATES catalog is as 
uniform as possible for the high-latitude (|6| > 10°) sky, limited by gaps in which the initial 
4.8 GHz data are unavailable. We believe that this catalog is an excellent starting point for 
comparison with other all-sky samples (e.g., microwave and 7-ray). 

Here we wish to find EGRET-like blazars, so we adopt the FoM of SRM03, which was 
derived from comparing the well-established 3EG blazar sources with the northern (CLASS- 
generated) subset of the CRATES catalog. This FoM is given by the heuristic fitting formula 
F0M3EG = 100 x P a x P s x P x x -Ptsj where the P terms are "excess probabilities" for the 
observed parameters for radio sources near 3EG sources. Here, P a = 0.19 — 0.35ai ow / 8 . 4 
(0 < P a < 0.4), P s = -3.47 + 2.451og 10 S 8 . 4 - 0.3 41og? n S* A (0 < P s < 1) and P x = 



0.99 + 0.411og 10 F (0.5 < P x < 1), with F the RASS flVoges et all 119991 ) counts per second 
and the P terms bounded to the ranges in parentheses. Finally, Pts = 1 — CL, where CL is 
the confidence limit of the 3EG source localization contour passing through the position of 
the radio source. In essence, the FoM is composed of the product of the "excess" probabilities 
of sources of a given flux density, spectral index, etc. over random chance. While the FoM 
probability is not directly normalized, "false positive" rates were computed at each FoM 
level by comparison with the statistics of scrambled versions of the sky catalogs. Of course, 
once we have an initial survey of LAT blazar sources, it will be appropriate to derive new 
coefficients, "re-training" the FoM against this sample. 

To develop an all-sky survey of blazar candidates, we compute an FoM for each source 
in the CRATES catalog. We must do this without reference to 3EG sources. Thus for this 
paper we define FoM = P a x Pg x Px- To connect with FoM 3 eg ; note that a blazar with 
the present FoM = 0.2 would correspond to a FoM 3 eg = 1 at the 95% localization contour 
of a 7-ray source, a "likely" (>90% correct) identification. With this definition, 5059 of the 
CRATES sources have a nonzero FoM. To focus our follow-up on the best and most 3EG-like 
objects, we define CGRaBS as those 1625 sources with FoM > 0.04. This corresponds to an 
SRM03 value of FoM 3EG = 2 at the 50% localization contour, a very likely association, and 
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a FoM 3 eg = 0.2 for a source at the 95% confidence contour, a reasonable (>80%) likelihood 
of an association. Figure 1 shows an Aitoff equal-area projection of the CGRaBS sample 
along with its parent survey, CRATES. Figure 2 shows a projection of the CGRaBS sample 
indicating the FoM of each source. 




Fig. 1. — Aitoff equal-area projection of the CRATES parent sample (small dots) and the 
CGRaBS sample (large dots) in Galactic coordinates (/, b). The central meridian is I = 0°. 
A few small holes are visible just below 5 = 0° (dot-dashed line), stemming from incomplete 
PMN sky coverage. 



The radio spectral index is a major component of our FoM; thus, since the interfero- 
metric observations at 8.4 GHz and low frequency were non-simultaneous, variability can, 
in principle, affect our FoM measur ements. Luck ily, the variability in the radio is modest 
compared to the high energy bands: iHealey et al.l found that the mean 8.4 GHz variability is 
<14% and the low- frequency variability on the relevant several- year timescale is even smaller. 
Thus, we do not expect that radio variability will dramatically affect our FoM estimates. 
Further, the (more likely variable) RASS detections turn out not to be a major selection 
bias in this survey. If the X-ray contribution to the FoM is ignored and a purely radio-based 
FoM is computed, then 98.5% of our sources still satisfy the CGRaBS FoM cutoff. Thus, 
while the X-ray flux from a small number of sources boosts them into the sample, the main 
effect of the X-ray contribution is to shuffle the ranking within the set of sources that are 
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Fig. 2. — Aitoff equal-area projection of the CGRaBS sample in Galactic coordinates (/, b). 
The central meridian is I — 0°. The radius of each dot is proportional to the FoM of the 
source; the dot for a source with FoM = 0.1 is shown for comparison. The dot styles indicate 
optical classifications (see §§3.2.1-2): • = FSRQ, <g> = BLL, © = AGN, o = unknown. 

already qualified. Since the radio FoM weighting increases for bright and inverted (rising) 
spectra, its net effect is to impose an effective extrapolated flux density limit at a higher radio 
frequency. For example, the FoM = 0.04 cutoff corresponds to an extrapolated flux density 
at 100 GHz of Sioo > 230 mJy (although we do not expect all sources to have a constant 
a to such high frequency). Less than 1% of the full CGRaBS targets have an extrapolated 
flux below this threshold, and these are all low- FoM sources with very high X-ray flux (i.e., 
largely high-peaked sources; see the next section). 

Three CGRaBS sources warrant special comment. The CRATES entry for J0352— 2514 
is a combination of 8.4 GHz observations at two epochs, one with an unflagged mapping error 
and a grossly erroneous position. The CGRaBS entry for J0352— 2514 uses only the good 
epoch to determine the correct position, the 8.4 GHz flux density, and the spectral index. 
Sources J0805— 0111 and J1639+1632 have nominal CGRaBS spectral indices (and thus 
FoMs) that are almost certainly overestimates. Their NVSS counterparts have marginally 
resolved jet structure, and the NVSS decompositions offset the core toward the jet. A 
faint, spurious counter-jet component was introduced and, being slightly closer to the 8.4 
GHz position, was selected as the 1.4 GHz counterpart, leading to a highly inverted spectral 
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index and a high FoM. We include these sources in the survey since they satisfy the CGRaBS 
prescription; a more careful treatment of the NVSS counterparts would give a smaller spectral 
index and FoM. This effect is quite rare, occurring in CGRaBS for only these 2 sources (out 
of 1625, or 0.12%) and in CRATES for no more than 20 sources (out of 11,131, or 0.18%). 



3. Optical follow-up 

We have specifically not required a previous optical (or X-ray) detection of our blazar 
candidates. This radio-driven selection allows us to sample completely the flat-spectrum 
sources and avoid biasing the detected population. For example, X-ray-bright sources are 
prefere ntially low-power "blue" b lazars such as BL Lacs (so-called high-peaked blazars, or 



HBLs; iPadovani fc Giommilll995l ). Similarly, requiring optically bright counterparts can bias 
the sample toward low redshift. However, since the principal goal of the CGRaBS project 
is to secure optical identifications, we do need good magnitude estimates. To maximize 
uniformity, we are working toward complete identification for R < 23. In practice, we have 
also observed a number of radio-bright and X-ray /7-ray-bright but optically faint sources 
beyond this limit to explore the extrema of the population. 



3.1. Counterparts and photometry 

One defining blazar characteristic is rapid optical variability. Thus, we must set a fiducial 



epo ch" for the optical magnitudes. In practice, we use the USNO-B1 catalog (IMonet et al. 



20031 ) since this is the largest source of suitable R magnitudes; we take the fiducial magnitude 
to be that of the more sensitive second epoch survey (R2). Since we have precise radio 
positions for the cores of all sources, we identify a USNO-B1 source as the counterpart of 
a CGRaBS source if the optical position is within 1.5" of the radio position. This gives a 
large fraction of the required magnitudes, with completeness dropping between R ~ 20 and 
R ~ 21. For the north G alactic cap, we can supplement these with SDSS identifications 



(through Data Release 5, lAdelman-McCarthy et al.l 120071 ) to r' < 22. In confused cases 



these archival data were examined visually to determine the best counterpart match. In a 
number of cases, we were also able to see clear counterparts that were too faint for inclusion in 
the USNO-B1 catalog but whose magnitudes could be reasonably estimated by measurement 
of the digitized plate data. In view of the variable blazar magnitudes and non-stellar colors, 
this low-precision photometry is adequate to guide the follow-up spectroscopy. 



To complete the process of optical identification (and to improve a few poor USNO-B1 
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magnitudes), we have conducted our own imaging campaign, primarily at the 5 m Hale 
Telescope at Palomar, the 3.6 m New Technology Telescope (NTT) at La Silla, and the 2.7 
m Harlan J. Smith Telescope at McDonald. Typical exposures were 180 s through Gunn 
r' under varying conditions, and magnitudes were calibrated against multiple field stars. 
For some particularly interesting sources (e.g., high radio-to-optical flux ratio candidates 
for high redshift), these were supplemented with izJHK imaging. We do not report here 
on these optical/IR SEDs. All r' magnitudes have been converted to R using the average 
color term (R — r' — —0.253) of CGRaBS sources detected by both the SDSS and USNO- 
Bl. A magnitude (or limit) for each source is listed in Table 2. For some of the lowest 
redshift sources, the magnitude is dominated by the flux from the (extended) host galaxy. 
We also list the nominal Gal a ctic e xtinction for the source direction Ar, derived from the 



Schlegel. Finkbeiner fc David (119981 ) maps. Even though the sources are at high latitude, 
there are a few targets behind dust clouds, indicating a large nominal extinction. However, 
we do not expect extinction to bias our measured population as the large Ar are not prefer- 
entially associated with the faint targets. Furthermore, only 4% of the blazars have Ar > 1 
and 0.5% have A R > 2; only 4 sources are excluded from the targeted R = 23 sample by the 
known extinction. As of 2007.5, there are 88 objects (5.4%) that do not have measured R 
magnitudes; of these, 45 have limits fainter than R = 23 and thus do not nominally require 
spectroscopy for the complete survey. The sources with brighter limits will be the subject 
of further imaging. Note that with 68 CGRaBS sources known to be fainter than R = 23, 
we expect that the survey will be >95% complete at this magnitude limit. 

Figure 3 shows the distribution of R magnitudes and limits. Since the Ar are in general 
small, the extinction-corrected histogram is very similar. At first sight, the rapid drop 
between 20 < R < 21 would seem to be due to the USNO-B1 survey completeness limit. 
However, we have sufficient deeper CCD imaging to determine that the drop in numbers is 
largely intrinsic, although we need to complete the imaging before we can characterize the 
details of the faint source distribution. The right panel shows that we need to complete 
identifications to faint magnitudes (R > 19) to get a representative sample of the higher 
redshift sources. 



3.2. Spectroscopy 

Our spectroscopic goals are a basic classification of the AGN type, redshift measurement, 
and measurement of emission line equivalent strengths and kinematic widths (for luminosity 
and mass evolution studies). Thus, the bulk of our new observations have been low- resolution 
1Z ~ 500—1500 long-slit spectroscopy. Most of the sources are flat-spectrum radio quasars 
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Fig. 3. — Left: Magnitude distribution of CGRaBS sources. Lower limits on R are shown 
by the dashed histogram. Right: Magnitude distributions for low-redshift and high-redshift 
sources. 

whose broad lines allow easy identification with relatively low signal-to-noise. However, a 
significant fraction of the sources (~15%) are weak-lined BL Lac sources. For these, we 
require high S/N and/or high resolution to determine the redshift from host absorption 
lines. Such measurements require long exposures with large telescopes. At present, we have 
identified sources as BL Lacs to R ~ 20, but our ability to measure the redshift drops 
significantly above R ~ 18.5; these BL Lacs are the subject of further spectroscopy at 
higher dispersion. In this paper, we present a progress report on the optical identifications. 
Additional papers will discuss the properties of the complete sample, the source SEDs, and 
the constraints on blazar evolution. 



3.2.1. Observations 



A fair fraction of the CGRaBS sources are bright, well-known AG N; thus, we have vetted 
our c atalog against the twelfth edition of t he Veron quasar catalo g (jVeron-Cetty fc Veron 
20061 ) and the SDSS DR5 quasar catalog (ISchneider et al.l 120071 ). We have also queried 
NED (http : //nedwww. ipac . caltech.edu/ ) for all CGRaBS sources to find any other red- 
shifts and identifications in the literature. Archival data identify ~45% of the CGRaBS 
objects (~60% of the redshifts in hand); the remainder are the targets of our own spec- 
troscopic campaigns. The great workhorse of our spectroscopic effort has been the 9.2 m 
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Hobby*Eberly Telescope (HET) at McDonald, which has observed hundreds of CGRaBS 
sources in the accessible declina tion band —11° < 5 < +73°. The telescope is fully queue- 



scheduled ( IShetrone et al.l 120071 ). allowing us to receive data remotely year-round and to 



spread the cost of inclement we ather and unfav orable conditions. We use the Marcario Low 



Resolution Spectrograph (LRS; lHill et al.lll998l ) with grism Gl (300 lines mm x ), 2" slit, and 



a Schott GG385 long-pass filter for a resolution of 1Z ~ 500. Typical exposures are 2 x 600 
s, providing redshifts of emission-line objects to R ~ 22; brighter objects are also observed 
under poor conditions with 2 x 300 s. 

In addition to our ongoing HET observations, we have mounted dedicated campaigns 
at a number of other facilities. We conducted three runs totaling 13 nights (over half lost 
to weather) on the 2.7 m Harlan J. Smith Telescope at McDonald, using the Imaging Grism 
Instrument (IGI) and the 6000 A VPH grism. We observed 28 objects with the 1.5 m 
telescope at Cerro Tololo in the "13/1" setup (grism 13, 150 lines mm -1 ) in service observing 
mode as part of the Small and Medium Aperture Research Telescope System (SMARTS) 
program. We conducted two runs totaling 8 nights on the 3.6 m NTT at La Silla with the 
ESO Multi-Mode Instrument (EMMI) in the low-resolution spectroscopy (RILD) mode and 
grism 2 (300 lines mm -1 ). To date, we have had three runs totaling 12 nights on the 5 m Hale 
Telescope at Palomar with the double spectrograph (DBSP), using a 300 lines mm -1 grating 
on the blue side and a 316 lines mm -1 grating on the red side. We observed 12 objects with 
the 8.2 m Kueyen telescope (the second unit telescope at the Very Large Telescope, VLT) in 
service observing mode with Focal Reducer/Low-Dispersion Spectrograph 1 (FORS1) and 
grism GRIS_300V (300 lines mm -1 ). Finally, we have had three runs totaling 4 nights on the 
10 m Keck I Telescope at Mauna Kea (however, the night of 2006 October 28 was the first 
observing night after the 2006 earthquake, and pointing was severely restricted; observations 
remained substantially constrained even on the night of 2006 November 24). For these 
observations, we used the Low Resolution Imaging Spectrometer (LRIS), employing a 600 
lines mm -1 grism on the blue side and a 300 lines mm -1 grating on the red side. A summary 
of the observations is shown in Table 1. 

The 1.5 m telescope observations were taken with a fixed N-S slit within a few hours of 
culmination. For all other systems, observations were taken with a long slit at the parallac- 
tic angle. Basic reduction steps were applied to the spectra using standard IRAF routines. 
Although every effort was made to minimize differential slit losses, in view of the variable 
slit widths and seeing, we have not attempted to derive absolute spectrophotometry. Af- 
ter standard star calibration, we estimate that the relative spectrophotometric accuracy is 
~30%, based on comparisons of observations of individual targets at different epochs with 
different instruments. Spectra were corrected for telluric absorption, and all observations 
for a given target were combined, weighted by S/N, to produce a final spectrum. Sample 
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spectra are shown in Figure 4. 




4000 5000 6000 7000 8000 4000 5000 6000 7000 8000 

Wavelength (A) Wavelength (A) 



Fig. 4. — Sample CGRaBS spectra. 



3.2.2. Results 



Our spectral analysis starts with a basic source classification. The vast majority (84%) 
are flat-spectrum radio quasars (FSRQs) dominated by strong broad emission lines. The 
weak-lined BL Lac (BLL) class is somewhat heuri stically defined; he re we designate as 
BLLs sources that exhibit the following properties (jMarcha et al.l 119961 ): (1) emission line 
equivalent width < 5 A, and (2) H/K 4000 A break contrast = (/+ - f~)/f + < 0.4, 
where / + (/~) is the flux density redward (blueward) of the break. It is often possible to 
establish that a source is a BLL even when the redshift is impossible to determine. For 
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Table 1. Summary of CGRaBS observations. 



Telescope 


Dates 


Wavelength 
range 

(A) 


Spectral 
resolution 

(A) 


Typical 
seeing 
(arcsec) 


Typical 
exposure 

(<0 


9.2 m Hobby*Eberly 


Ongoing, 2002— present 


4100-9700 


17 


1.5 


600, 1200 


2.7 m Harlan J. Smith 


2005 May 27-31 
2005 Oct 27-31 
2007 Mar 26-28 


4250-8250 
4250-8250 
4250-8250 


12 
12 
12 


1.5 
2.0 
3.0 


600, 1200 
600, 1200, 1800 
600, 1200, 1800 


1.5 m CTIO 


2005B 


3500-9000 


17 


1.5 


1200, 1800 


3.6 m NTT 


2006 Aug 29-Scp 1 
2007 Jan 22-25 


3900-9100 
3900-9100 


10 
10 


1.3 
1.0 


600, 1800 
600, 1200 


5 m Hale 


2005 Nov 5-9 

2006 Aug 17-18 

2007 Jan 15-16 
2007 Apr 19-21 


3300-9500 
3300-9500 
3300-9500 
3300-9500 


5 a , 16 b 
5 a , 16 b 
5 a , 16 b 
5 a , 16 b 


2.5 
1.7 
2.0 
2.5 


600, 1200 
600, 1200 
600, 1200 
600, 1200 


8.2 m VLT-Kueyen 


Period 78 


3500-8000 


17 


1.2 


600, 900, 1200, 1800 


10 m Keck I 


2006 Jul 22-23 
2006 Oct 28 
2006 Nov 24 


3300-9300 
3300-9300 
3600-9600 


3 a , ll b 
3 a , ll b 
3 a , ll b 


1.5 
1.5 
2.5 


600, 1200 
600, 1200 
600, 1200 



a Blue side. 
b Red side. 
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sources with R > 15, we compute M R for the ACDM concordance cosmology (smaller R 
values are usually host-dominated, in any case) and classify broad emission line sources with 
Mr > —23 as AGN. Thus, we list here three blazar designations: continuum-dominated 
BLLs, low-luminosity broad-line AGN, and luminous broad-line FSRQs. A small number of 
non-blazar sources is also present. Sources with narrow lines (f fwhm < 1000 km s _1 ) are 
denoted as narrow-line radio galaxies (NLRGs). Sources with small line equivalent widths 
but large H/K break contrasts are denoted as galaxies. These low-redshift sources may 
represent the low-luminosity extension of the blazar phenomenon. One extremely compact 
planetary nebula (PN) made our survey cuts. Finally, in four cases, the radio position was 
within 1.5" of a field star that dominated the initial spectrum. With improved imaging, the 
fainter CGRaBS blazar counterparts can be identified. 

Redshifts were measured by cross-correlation analysis. For a modest number of FSRQs, 
only a single strong, broad emission line is identified. In most cases, we conservatively identify 
this as Mg II A2800, supported by the absence of strong lines expected for other identifications 
and, often, by Fe II structure in the surrounding continuum. Nevertheless, these redshifts are 
flagged by a colon (:), indicating possible systematic uncertainty. Absorption line redshifts 
were obtained for some BLLs. In a few cases, the BLL sources had multiple observations, 
and we were able to obtain emission line redshifts when the source was in a low continuum 
state and the emission line equivalent widths were relatively large. A few additional BLLs 
have redshift constraints, with upper limits from the lack of Ly-a absorption in the UV and 
lower limits from clearly identified (typically Mg n) intergalactic absorption systems. We 
have also measured continuum flux densities and equivalent and kinematic widths for the 
strong optical/UV resonance lines. These will be used to study the black hole masses and 
evolution. 

Table 2 presents the first page of the CGRaBS catalog; the full table appears in the 
online edition. Here we include the precise position, the 8.4 GHz core flux density, the FoM, 
the R magnitude, the extinction Ar, and the optical classification and redshift, if any. 

4. Discussion 

To date, we have 1226 redshifts and 64 BLLs with unknown redshift. Thus classification 
is 79% complete with respect to the entire survey and 85% for objects with known R < 23. 
So far, 10.3% of all objects classified are BLLs, 3.4% are AGN, and 1% are NLRGs. Figure 
5 shows the completeness as a function of magnitude. Source classification and redshifts are 
>85% complete to R = 20. While the completeness drops off rapidly beyond this, so do 
the source counts, and so reaching >95% completeness at the survey limit is feasible. Note, 
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Table 2. The CGRaBS catalog. 









58.4 GHz 


X-ray R 


Ar 






Name 


R.A. a 


Decl. a 


(mjy) 


FoM fiag b (mag) 


(mag) 


Type c 


z 



J0001 — 1551 


00 


01 


05 


.3.3 


— 15 


51 


07. 


1 


335 


.90 





050 




18 


09 





.08 


FSRQ 


2 


.044 


J0001+1914 


00 


01 


08 


.62 


+19 


14 


3.3 


8 


504. 


.20 


0. 


,105 




20 


,50 





,11 


FSRQ 


3 


.100 


J0003+2129 


00 


03 


19 


35 


+21 


29 


44. 


.4 


269 


.70 





.098 




19 


75 





12 


AGN 





.450 


JU004— 1148 


00 


04 


04 


.92 


— 11 


48 


58. 


.4 


771. 


.90 





112 




19 


09 





08 


BLL 






J0004+4615 


00 


04 


16 


1.3 


+46 


15 


18. 





214 


.80 





060 




20. 


.44 





,24 


FSRQ 


1 


.810 


J0004+2019 


00 


04 


35 


76 


+20 


19 


12 


2 


162 


.50 





.058 




20. 


.25 





,10 


BLL 





.677 


J 0004— 4736 


00 


04 


35 


.68 


—47 


36 


18 


6 


780 


.40 





.067 




15 


.88 





05 


FSRQ 





.880 


J 0005— 1648 


00 


05 


17 


.9.3 


— 16 


48 


04. 


7 


281 


.70 





.050 




18 


.37 





,07 








J 0005+0524 


00 


05 


20 


.21 


+05 


24 


10. 


.7 


228 


,90 


0. 


.066 




16. 


.26 





,08 


FSRQ 


1 


.900 


J0005+3820 


00 
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b "X" indicates that a source would not satisfy the FoM cutoff if its X-ray flux were ignored. See §2. 
c See §§3.2.1-2 for discussion of the type classifications. 
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however, that only ~52% of the BLLs have redshifts and that this fraction falls off quickly 
above R = 18. Clearly, pushing the largely complete BLL sample fainter than R = 20 will 
be a challenge. 
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Fig. 5. — Completeness as a function of R for source identification and redshift measurement 
(histograms). The points show the fraction of identified BLL with measured redshift; small 
numbers lead to substantial error bars. 

We defer full discussion of the sample properties until we reach our expected 95% com- 
pleteness to R = 23. However, it is already interesting to examine the redshift distribution 
of the sources detected to date (Figure 6). The non-BLL (largely FSRQ) distribution peaks 
at z ~ 1.3 and has an exponential fall-off (dN/dz oc 10~ a6z ) to high redshift, extending to 
z = 5.5. From SED information on optically faint sources, we expect the high-redshift pop- 
ulation to increase somewhat in the complete CGRaBS sample, but it is clear that there will 
be only a handful of radio-bright blazars at z > 4. If any of these are detected by the LAT, 
as expected, they will be particularly important targets for multiwavelength spectral and 
variability studies. In fact, with only ~40 sources at z > 3, careful study of these few high- 
redshift objects will be important for several LAT programs, e.g., extragalactic background 
light (EBL) studies and studies of jet evolution and interaction with the CMBR. 
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Fig. 6. — Redshift distributions for the (partly complete) CGRaBS survey. The solid-line 
histogram shows FSRQs. The short-dashed histogram gives the redshift distribution for 
solved BLLs. The long-dashed histogram shows a variety of other AGN (NLRGs, passive 
ellipticals, etc.), which contribute only at very low redshift. 



We are also assembling an important new sample of radio-bright BLLs. To date, we 
have 133 sources definitively classified as BLLs, but this will likely grow since a substantial 
number of other sources have observed BLL-like spectra but need somewhat improved S/N 
observations to exclude emission lines with EW > 5 A throughout the observed spectrum. 
Among the brighter sources ~15% are BLL; at this incidence, we expect ~245 sources to 
have a final BLL classification. As noted, it will be very tough to obtain redshifts of the 
faintest BLLs. However, the 70 redshifts already in hand rep resent a substantia l radio-bright 
sample. For example, it is twice the size of the 1 Jy sample (jStickel et al.lll99ll ) and extends 
to nearly twice the redshift. At present, we have 11 BLLs at z > 1, a third of all known 
z > 1 BLLs, so that the full survey should be useful for probing evolution of this population. 



Of course, the most important application of the CGRaBS catalog is the identification 
with other all-sky samples and the generation of multiwavelength SEDs. We are already ex- 
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amining the radio to X-ray spectra of these sources and eagerly look forward to the upcoming 
sky surveys with AGILE, the air-Cerenkov TeV observatories, and especially GLAST, which 
will measure the 7-ray power peak expected for many of these sources. 
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